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ABSTRACT 

The M82 'cap' is a gas cloud at a projected radius of 11.6 kpc along the minor axis of this well 
known superwind source. The cap has been detected in optical line emission and X-ray emission and 
therefore provides an important probe of the wind energetics. In order to investigate the ionization 
source of the cap, we observed it with the Kyoto3DII Fabry-Perot instrument mounted on the Subaru 
Telescope. Deep continuum, Ha, [N II]A6583/Ha, and [S II]AA6716,6731/Ha maps were obtained with 
sub-arcsecond resolution. The superior spatial resolution compared to earlier studies reveals a number 
of bright Ha emitting clouds within the cap. The emission line widths (< 100 km s~ x FWHM) and 
line ratios in the newly identified knots are most reasonably explained by slow to moderate shocks 
velocities (u s hock = 40-80 km s _1 ) driven by a fast wind into dense clouds. The momentum input 
from the M82 nuclear starburst region is enough to produce the observed shock. Consequently, earlier 
claims of photoionization by the central starburst are ruled out because they cannot explain the 
observed fluxes of the densest knots unless the UV escape fraction is very high (/ esc > 60%), i.e., an 
order of magnitude higher than observed in dwarf galaxies to date. Using these results, we discuss the 
evolutionary history of the M82 superwind. Future UV/X-ray surveys are expected to confirm that 
the temperature of the gas is consistent with our moderate shock model. 

Subject headings: galaxies: individual (M82) — intergalactic medium — galaxies: ISM — galaxies: 
starburst 



1. INTRODUCTION 

Superwinds are galaxy scale outflows, caused by su- 
pernovae in nuclear starburst regions or active galac- 
tic nuclei (AGNs). They are so powerful that interstel- 
lar matter within the galaxies is blown out. Some of 
the material may escape to the intergalactic or group 
medium, while some of the material may be re cycled 
throughout the galactic halo ([Cooper et al.l 12008). Su- 
perwinds are expected to quench star-formation activity 
(feedback) and to enrich the external medium with new 
metals. Generally, galactic winds are diffuse and difficult 
to observe. M82, one of the neare st starburst galaxies 
(3.63 Mpc, iFreedman etaLl fl99l . is one of the most 
well known examples of the sup erwind phenomenon. Its 
large inclination ang le (i ~ 80°: iLynds fc Sandagdl 19631 : 
iMcKeith et al.lll995j l and proximity allow us to see many 
details of the wind phenomenon far from the galactic 
plane. The source has been observed in hot gas (T ~ 10 6 
K; e.g.. lCappi et al.lll999f ). ionized gas (T - 10 4 K; e.g., 
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IMcCarthv et al.lll987t iBland fc Tullvi[l988l), and molec- 
ular g as (T ~ 10 2 K; e.g.. lNakai et al.lll987tlWalter et al.1 
120021 . 

The kinematics and ionization of the wind ma- 
terial over the inner few kiloparsecs h a ve be en in- 
vestigated in detail. IMcKeith et al.l (|T995) and 
iShopbell fc Bland-Hawthorn! (|1998l ) modeled the out- 
flow structure using position-velocity diagrams in op- 
tical emission lines. The emission line ratios of the 
inner region indicate that photoionization by the nu- 
clear starburst plays a significant role in the excita- 
tion ^McCarthy et aljl987t IShopbell fc Bland-Hawthorn] 
1998). In recent years, new observ ational methods such 
as inte gral field spectroscopy (e.g.. iWestmoauette et al. 
2009al lbr and spectropolarimetry (e.e.. lYoshida et~ 
201 lh have revealed its more intricate structure. 

Our goal is to shed light on processes behind large-scale 
galactic winds. Very little is known about their total ex- 
tent, energetics and importance in the context of galaxy 
evolution. By studying the most spatially extended emis- 
sion, we can obtain a better understanding of the total 
kinetic energy of the wind. There are many questions 
that remain unanswered for M82's outflow. How old is 
the wind and how far does it extend? Is it powered by 
radiation pressure or wind pressure, or a combination of 
both? Is the source of energy impulsive or sustained over 
many dynamical times? Is most of the outflowing mate- 
rial swept up or entrained from the disk? Does the wind 
material escape the galaxy or fall back to the disk? To 
have any chance of answering these questions, we need a 
better understanding of the most basic properties of the 
large-scale wind. 

The most distant gas cloud in M82 is the 'cap' origi- 
nally discovered in Ha and X-ray emissi on at a radius of 
11.6 kpc along the minor axis of M82 CD cvmc fc Ballvi 
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Il999t iLehnert et al.lll999l) . Strong UV emission provides 
evidence fo r reflecting dust in the cloudlets that make 
up the cap (|Hoopes et al.ll200"5l) . The metal abundances 
of O, Ne, Mg, Si, and Fe of X-ray emitting gas in the 
cap suggest that most of the metals arise from a cir- 
cumnuclear starbu rst dominated by Type II supernovae 
(jTsuru et all 12051 . 

We now show that the dominant ionization source in 
the cap provide s an importan t clue to the wind's ori- 
gin and history ILehnert et all (|1999| ) suggested the cap 
is either photoionized by UV photons from the nuclear 
starburst region or by a shock being driven by the hot 
wind into a dense halo cloud, or a combination of both. 
The X-ray observations already support the idea that the 
wind reaches the distance of the cap, but are the optical 
emission line diagnostics consistent with a wind-driven 
shock? 

Therefore, in order to obtain emission line intensity 
map and line ratio maps at high spatial resolution, we 
carried out Fabry- Perot observations of M82's cap with 
the Subaru Telescope. This combination enables us to 
detect weak emission with a larger field of view than 
that of integral field spectroscopy. Through a compar- 
ison of the observed line ratios and those calculated by 
photoionization and shock models, we discuss the ioniza- 
tion source of the M82 cap and a likely evolution history 
for the large-scale galactic wind. 

2. OBSERVATIONS AND DATA REDUCTION 

On 2011 November 22, we observed the central part of 
the M82 cap, -10' N of the nucleus of M82, wi th the Ky- 
oto3DII Fabry- Perot mode (jSugai et all 120101 ) mounted 
on the Cassegrain focus of the Subaru Telescope. Fig- 
ure [1] displays the position of the cap relative to the 
M82 center, and indicates the region where we observed 
in this observation. This mode uses an ET-50 etalon 
manufactured by Queensgate Instruments. The field 
of view is ~1'.9 and the pixel scale is 0".112 pixel -1 
after 2x2 on-chip binning. The spectral resolution 
R 348 corresponds to 19 A at 6598. 95A. We obtained 
14 object frames for Ha + [N II]AA6548,6583, five for 
[S II]AA6716,6731, and two for the off bands. 

The observed wavelengths at the field centres are sum- 
marized in Table [TJ The exposure time for each frame 
was 300 seconds. We also o bserved a standard sta r 
EGGR247 for flux calibration (jBland-Hawthornl IT995I) . 
Bias subtraction and flat fielding were performed for the 
target and standard star frames. Because the center 
wavelength in Fabry-Perot observations depends on the 
distance from the center of the field of view, simple sky 
subtraction results in some residuals due to sky emission 
lines. We measured sky emission fluxes in blank regions 
of the object frames, and subtracted it from the regions 
at the same distance from the center. Flux calibration 
and distortion correction were carried out for t h e tar - 
get frames. We used a spectrum catalog of lOkd (1990) 
for flux calibration for each wavelength setting. The po- 
sitional offsets among the object frames were detected, 
because the Cassegrain auto guider was unavailable due 
to repairs and we did not use it in this observation run. 
We corrected the offsets by using the stars in the target 
frames. We matched the spatial resolution of the tar- 
get frames to the worst one, 0".9, and carried out 4x4 
binning, resulting in the pixel scale of 0".45 pixel -1 . 
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Figure 1. Ha (greyscale) and X-ray (contours) images of the 
whole M82 ULehnert et ail 119991 ). The field of view of Figure [2] 
panels is displayed with a rectangle. 



Table 1 

Observed wavelengths for the M82 cap 



band 



wavelength [A] 



Ha + [N II]AA6548,6583 
[S II]AA6716,6731 
continuum 



6546, 6554 a , 6562 a , 6570 a , 6578 a , 6586 a , 6594 a 
6714, 6722, 6730, 6738, 6746 
6656 a 



661 



Note. — These are the observed wavelengths at the center of our field of viei 
a Two exposures were performed at these wavelengths. 

Because of the relatively low spectral resolution, Ha 
and [N II]AA6548,6583 were blended in Ha+[N II] band. 
We fitted these lines pixel by pixel with the transmis- 
sion c urve of the Fabry-Per ot interferometer (Airy func- 
tion: iBland fe Tullylll989l) and decomposed them (see 
iMatsubavashi et all 12009). Better wavelength sampling 
than that of the previous observation enables us to find 
the best velocity center. We fitted the emission line fluxes 
at each velocity, from 100 km s -1 to 700 km s -1 , and se- 
lected the velocity for which the fitting residual is the 
smallest. For line decomposition, we assumed that the 
[N II]A6548 flux is one-third of [N II]A6583 flux, and 
that the velocity centers of Ha and [N II]AA6548,6583 
are same. 

The velocity dispersion is fixed to km s -1 because 
it is much smaller than the spectral resolution of the 
instrument. This assumption is reasonable, since the 
observed veloci ty dispersion of Ha a t the cap is small 
(~100 km s -1 ; iDevine fe BaUvl ll999). The same fitting 
was performed for the [S II] band data. However, the 
[S II]A6716/[S II]A6731 cannot be determined well, be- 
cause the wavelength difference between these two lines 
is smaller than the spectral resolution in this observation. 
Therefore we only use the total flux of [S II]AA6716,6731 
in this study. Our la detection limit in Ha surface 
brightness is estimated to be 6.5 x 10 -18 erg cm -2 s -1 
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Table 2 

flux and line ratios at each knot 



Figure 2. (a) Continuum surface brightness, (b) He* intensity, (c) 
[N II]A6583 intensity, (d) [S II] AA6716,6731 intensity, (e) [N II] /Ha 
ratio, and (f) [S II] /Ho ratio maps of a part of the M82 cap. The 
color scales are 0-3.3 X 10 — 19 erg cm -2 s -1 A -1 (0.45 arcsec) -2 , 
0-2.5 X 10 -17 erg cm -2 s -1 (0.45 arcsec)" 2 , 0-1.5 X 10 -17 erg 
cm -2 s -1 (0.45 arcsec) -2 , 0-1.5 X 10 -17 erg cm -2 s -1 (0.45 
arcsec) -2 , 0-1.0, and 0-1.0 for the continuum, Ho, [N II], [S II], 
[N II] /Ho, and [S II] /Ho maps, respectively. Two bright stars ap- 
peared at upper side in the panel (a) are masked in the panels 
(b)-(f). North is up and east is left. Bar in lower left in each 
panel represents 30", which corresponds to 0.5 kpc. Knot IDs arc 
displayed in the Ho intensity map. (A color version of this figure 
is available in the online journal.) 



Knot ID 


Ho flux 


[N Il]A6583/Ho 


[S Il]AA6716,6731/Ho 




(10 -16 erg cm -2 s -1 ) 


E-C 


4.3 ± 0.1 


0.34 ± 0.03 


0.63 ± 0.06 


E-N 


4.0 ± 0.1 


0.22 ± 0.03 


0.60 ± 0.06 


E-El 


3.3 ± 0.1 


0.32 ± 0.04 


0.63 ± 0.08 


E-E2 


2.8 ± 0.1 


0.27 ± 0.04 


0.45 ± 0.08 


I - v. 


o.z in u.r 


u.io in u.uo 


u.ou in u.u/ 


F-N 


2.6 ± 0.1 


0.20 ± 0.04 


0.60 ± 0.09 


F-S 


2.5 ± 0.1 


0.20 ± 0.05 


0.55 ± 0.08 


G-C 


2.3 ± 0.1 


0.21 ± 0.05 


0.56 ± 0.10 


G-SEl 


3.1 ± 0.1 


< 0.11 


0.63 ± 0.08 


G-SE2 


3.1 ± 0.1 


0.12 ± 0.04 


0.62 ± 0.08 


G-SW 


2.9 ± 0.1 


0.18 ± 0.04 


0.52 ± 0.08 


H-C 


5.7 ± 0.1 


< 0.06 


0.55 ± 0.04 


H-W 


3.3 ± 0.1 


< 0.10 


0.66 ± 0.07 


J-C 


2.8 ± 0.1 


0.24 ± 0.04 


0.55 ± 0.09 


J-E 


2.9 ± 0.1 


< 0.12 


0.46 ± 0.07 


J-S 


2.6 ± 0.1 


0.17 ± 0.04 


0.50 ± 0.08 


Note. 


— The aperture sizes are 


2". 25 x 2". 25. 


Uncertainties are in 1 a 



levels, while upper limits are in 3 cr. 

is 5"— 10", which corresponds to 90-180 pc at the distance 
of M82. 

We renamed the Ha knots after iDevine fe Balfyl 
(1999), as shown in the Ha image (Figure [5] (b)). The 
Ha flux and luminosity of the brightest knot H-C in 
5" x 5" aperture are 1.7 x 10~ 15 erg cm" 2 s" 1 and 
2.8 x 10 36 erg s _1 , respectively. The electron density 
and mass of the brightest knot H-C are estimated as 1.0 



(3.3 cm b pc) at a temperat ure of 10 K. We ad opt a 
distance of 3.63 Mpc to M82 (jFreedman et al.lll994D . 

3. RESULTS 

Figure [5] (a) displays the continuum surface bright- 
ness map of the central part of the M82 cap. Only 
stars in the Galaxy and distant galaxies are detected. 
The number counts are consistent with the freely avail- 
able GalaxyCount program which provides source statis- 
tics for any window function down to 28 AB mag 
(jEllis fe Bland-Hawthornl 120071) . The relative posi- 
tions of objects in o ur image coincid e with those of 
Devine fe BaEvl (fl999T ) and SDSS DR7 (|Abazaiian et al.l 



arcsec 2 or an emission measure of roughly 1 Rayleigh fnl/i cm 3 ano - ^ x ^ 3 /Hai-^0i respectively, where 



-,3 f 1/2 



2009). In contrast, we cannot detect continuum emission 



from the cap. The upper limit in surface brightness at 
6656 A is 23.7 mag arcsec" 2 (AB, 5a), which corresponds 
to stellar mass of approximately 3 x 10 7 Mq using the 
cap size as 0.5 kpc 2 and mass-to-luminosity ratio at so- 
lar metallicity, star- formation history of SSP, and an age 
of 1 Gyr estimated from the iBruzual fe Chariot! (|2003f > 
model. This fact indicates that the cap is not a dwarf 
galaxy (jLehnert et al.| [l999'). 

In the Ha intensity map (Figure [2] (b)), clumpy 
and filamentary structures in the cap are clearly de- 
tected. Our Ha map is roughly consistent with that of 
IDevine fe Ballvl (| 1999(1 . but bright knots C and D identi- 
fied by them are not confirmed in our Ha image; instead 
these are detected in our con tinuum image. The Ha im- 
age of IDevine fe Ballvl ()1999D is clearly contaminated by 
continuum emission. Knots C and D are not related to 
the cap, and appear to be more distant disk galaxies. 
Our high resolution (= 0".9) Ha map enables us to re- 
solve some bright Ha knots. The typical size of Ha knots 



/hq,i indicates the filling factor of the knot. They are 
estimated from its Ha luminosity, size (5" = 90 pc), 
and the Ha recombination rate a ^a = 8.7 x 10 -14 cm 3 
s _1 ([Osterbrock fe FerT and 2006), with the assumption 
of spherical symmetry and a completely ionized gas. The 
observed Ha flux at each knot is displayed in Table [21 
The total Ha flux of the cap region in our field of view is 
7.3 x 10" 14 e rg cm" 2 s" 1 . Th i s flux is abo ut half of that 
estimated by IDevine fe Ballvl ()1999f l and iLehnert et al.l 
(|1999(1 (~1.3 x 10~ 13 erg cm -2 s -1 ) consistent with the 
fact that about half of the cap region falls within our 
field of view. 

Figures [1(c), (d), (e), and (f) displays [N II]A6583 in- 
tensity, [S II]AA6716,6731 intensity, [N II] /Ha ratio, and 
[S II] /Ha ratio maps of a part of the M82 cap, respec- 
tively. The [S II] flux map is generally similar to the 
Ha map. All counterparts of Ha knots are also found in 
the [S II] map. [S II] /Ha flux ratios are almost constant 
among these knots, 0.45-0.66. The [N II] flux map is 
quite different from the Ha flux map. Knots E-C and 
E-E are clearly detected in [N II], but the counterparts 
of the other knots, even the Ha brightest knot H-C, are 
barely detected. The observed [N II] /Ha flux ratios of 
knots E-C and E-El are the largest and peak at ~0.33. 
[N II] /Ha and [S II] /Ha line ratios at these knots are 
summarized in Table [2] and plotted in Figure [3] Small 
[N II] /Ha and [S II] /Ha ratios are consistent with the 
prev ious result of non-dete ction of forbidden lines at the 
cap ([Devine fe Ballvlll999l ). We do not find correlations 
between the Ha flux, [N II] /Ha and [S II] /Ha ratios. 

We compare the [N II] /Ha and [S II] /Ha ratios of 
the cap with those at other regions in M82. The emis- 
sion line ratios at a radius of 1 kpc from M82's center 
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[NII]/Ha 

Figure 3. Observed [S irj/Ha against [N Il]/Ha ratios listed 
in Table [2] Filled circles, open circles, and open triangles rep- 
resent the observed ratios at M82 cap (this study), superwind re- 
gions ~1 kpc from t he nucleus jShopbcll & Bland-Hawthor nlll998l; 
l Yoshida ct al. 201 1]), an d central regions l lO'Connell & M angano 
119781 ; ISmith et al.H2006D . respectively. Uncertainties in this study 
are in 1 a levels, while upper limits are in 3 a. Errorbars on the 
superwind region data represent the ranges of line ratios in their 
field of vie w. Line ratio grids calculated by a photoionization model 
(Cloudy: Fcrland ct al. 1998|) are plotted in solid lines. The as- 
sumed electron densities are (a) 1 cm -3 and (b) 2000 cm . The 
positive and negative values on the model lines denote metallic- 
ity and ionization parameter, respectively. Dashed line in panel 
(a) displays the line ratios calculated by a slow shock model of 
Shull & McKcc ( 1979) with solar mctallicity as a function of shock 
velocity. Data points for a shock velocity of 60, 80, and 90 km s — 1 
appear in this panel. The dotted lines in panels (b) and (c) dis- 
play t he line ratios calculated by a fast shock model of [Allen et all 
(2008) with solar and LMC metallicities, respectively, as a function 
of shock velocity at various magnetic field strengths. The left edge 
of each line corresponds to the lowest shock velocity, 200 km s — 1 . 

(|Shopbell fc Bland-HawthOT3l99cl lYoshida et al.ll20TTI) 
and in the circumnuclea r regions (jO'Connell fc Mangano! 
ll978tlSniith et al.ll2006ft . are also plotted in FigureEl We 
find an interesting trend of the line ratios with galactic 
radius. The [N II]/Ha ratio tends to decrease, while the 
[S II] /Ha ratio tends to increase with distance from the 
M82 nucleus. This fact suggests that some parameters, 
such as metallicity and shock velocity, gradually change, 
if the ionization source is the same in these regions. In 
FigureEl the line ratios of the M82 starburst regions and 
the cap are significantly different. Therefore, dust reflec- 
tion of the M 82 starburst region s, suggested by strong 
UV emission (Ho opes et al.l 120051) . is not the dominant 
emission mechanism in optical wavelength. 

Figure H] compares the observed [N II] /Ha and 
[S II] /Ha ratios of t he M82 cap with those of H II regions 
dJansen et alJl2000D. blue comp act galaxies (|Kong et al.l 
12002ft . LINERs (|Ho et al.lll997l) . a nd very extended ion - 
ized gas (EIG) in the Coma cluster (|Yoshida et al.|[2012T ). 
We find that the line ratios of the M82 cap are similar 
to those of some EIG knots which have larger [S II] /Ha 
ratios than the main sequence of star-forming galaxies. 
This fact suggests that the emission line ratios of the M82 
cap are not peculiar, and that the M82 cap and EIGs in 
the Coma cluster are ionized by the same mechanism. 

4. DISCUSSION 

First we discuss the dominant ionization source of 
the M82 cap. The observed [N II] /Ha and [S II] /Ha 
line rati os fall within the range of star-forming galax- 
ies (e.g..lLiang et al.ll2006|) an d diffuse ionized gas (e.g., 
iMoustakas fc Kenni cutt 2006]). The ionization source of 
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Figure 4. Observed [N II] /Ha vs [S II] / Ha of M82 cap (fille d 
circle) and EIGs in the Coma cluster (cross) ( Yoshida ct al. 2012). 
The shaded region re presents the line r atios of H II regions of 
nearby field gala xies (lanscn ct al. 2000), blue compact galaxies 
IIKong et ali r2002). and LINERs (Ho ct al. 1997). 



star-f orming galaxies is U V photons from massive stars 
(e.g., iKewlev et al.l l200lh . whereas there are some pos- 
sible sources for diffuse ion ized gas and EIG , such as 
photoionization, shock (e. g.. lAllen et al.ll2008D . and tur- 
bulent mixing layers fe.g., ISlavin et al.lll993D . In the case 
of turbulent mixing layers, the Ha emitting ionized gas 
should exist at the boundary of the hot gas, but they 
appear to be spatially coincident; thus, we rule out this 
model. In order to reveal the ionization source of the 
cap, we compare the observed line ratios with theoreti- 
cal values for photoionization and shock models. 

4.1. Photoionization 

Photoionization of the cap by the M82 nucle a r star - 
burst region was suggested by ILehnert et al.l ((1999). 
They calculated the number of ionizing photons at the 
cap region, and found that there are enough ionizing pho- 
tons relative to that estimated from Ha luminosity. This 
model assumes that ionizing photons are not strongly 
absorbed or scattered by interstellar matter between the 
M82 starburst regions and the cap region, presumably 
because the wind phenomenon has cleared the sight line 
of obscuring matter. 

The observed [N II] /Ha and [S II] /Ha r atios (Table 
H) ar e similar to those of H II galaxies (e.g. iLiang et al.l 
2006), and the observed ratios do not exceed the maxi- 



mum starburst line defined by IKewlev et al.l (|2001l) , un- 
less the [O III]A5007/H/3 ratio is larger than ~ 3. But 
the observed [S II] /Ha ratio is roughly equal to the 
largest value of ~ 40,000 SDSS star-forming galaxies 
(|Liang et a l. 2006) making UV photoionization less prob- 
able. 

In order to clarify whether photoionization can pro- 
duce the observed larger [S II] /Ha ratios, we calcu- 
lated e mission line ratios w ith a Cloudy photoionization 
model (jFerland et al.l [1998ft . We used two electron den- 
sities for the calculation: 1 cm" 3 for the cap, and 2000 
cm" 3 for the M8 2 center (jO'Connell fc Mangand 119781: 
iSmith et al.l 120061) . Due to its secondary nature, we as- 
sume that the nitrogen abundance scales with metallicity 
(Zg as ) as Zg as and that the stellar metallicity (Z*) is the 
same as Z gas (see iNagao et alj[201ll). T he SED of the 
stars is taken from iLeitherer et al.l (11999ft . 

Figure |3] (a) indicates that most of the observed line 
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ratios at the M82 cap are reproduced by the photoion- 
ization model of ionization parameter U = 10~ 35 -10 -30 
and metallicity Z ~ 0.4 Z@. U represents the dimen- 
sionless ratio of the ionizing photon density to the elec- 
tron density. To constrain t he ionization p a ramet er U of 
the M82 cap, we follow the iLehnert et al.l (|1999( ) condi- 
tion: the dimension of 3.7 x 3.7 kpc 2 , electron density of 
4.3 x 10 -2 /t^ 2 cm ~ 3 i where /h q ,2 indicates the filling 
factor of ionized gas averaged over the whole M82 cap, 
and the number of ionizing photons at the cap of 2 x 
10 50 photon s _1 . The calculated ionization parameter 

o 1/2 

U = 1.2 x 10 fa a 2 is consistent with that estimated 
from emission line ratios, if we assume /h q ,2 = 1- The 
larger distance from the ionization source than for typ- 
ical H II regions leads to smaller ionization parameter 
and larger [S II] /Ha ratios. 

The number of the ionizing photons which each knot 
receives is almost the same, and the ionization parame- 
ter depends on the knot electron density only. The cal- 
culated line ratios whose electron densities are smaller 
than 10 cm~ 3 are almost the same. For these reasons, 
we can consider that the difference in line ratios among 
the knots in the cap is explained by the difference in the 
electron density and / or metallicity, as seen in Figure |3] 
(a) . Although the lower density model cannot reproduce 
the ratios at the central regions, log U — —2.5 and Z = 
2.0 of the higher density model fits the observed ratios 
(Figure [3j (b)). 

However, taking account of the escape fraction of ioniz- 
ing photons, we find that this picture cannot be correct. 
ILehnert et al.l (|1999l ) uses Ha surface brightness aver- 
aged over the projected length of the M82 cap for the 
calculation of the escape fraction of ionizing photons. 
Since the required escape fraction is only ~3 %, they 
considered that the M82 starburst regions can provide 
enough ionizing photons for the M82 cap. Meanwhile, 
our high resolution image reveals that M82 cap is patchy 
and its filling factor /h q ,2 is much smaller than 1 (Fig- 
ure [2] (b)). This means that the Ha surface brightness of 
the brightest knot H-C must be explained by the ionizing 
photon flux from the M82 starburst region. Knot H-C is 
smaller and denser than the whole M82 cap, and there- 
fore much larger escape fraction is required. We estimate 
the required escape fraction of ionizing photons from the 
M82 starburst regions. We assume that the knots are 
spherical symmetric (Figure El). Then we can use the 
same calcula tion method asfBland-H awthorn fc Malonevl 
(19991 l200l . The Ha surface brightness of knot H-C is 1 
x 10~ 16 erg cm -2 s _1 arcsec -2 in 2". 25 x 2". 25 aperture 
(Tabled}. 

Since the surface brightness of 6 x 10~ 18 erg cm~ 2 
s _1 arcsec -2 corresponds to 1 Rayleigh at the Ha wave- 
length, the Ha emission measure of the knot H-C is 
16.5 Rayleigh. The required number of ionizing pho- 
tons is 3.7 x 10 7 cm~ 2 s _1 at the distance of M82, us- 
ing t he Ha recombination rate a na = 8.7 x 10 -14 cm 3 
s" 1 (jOsterbrock fc Ferlandl 120061) . Whereas the ioniz- 
ing photon flu x from the nuclear s tarburst region is 10 54 
photons s _1 (jMcLeod et al.l 119931 ). and 6 x 10 7 cm" 2 
s _1 at the distance of knot H-C from the M82 nuclear 
starburst region. Thus the required escape fraction for 
knot H-C is 60 %. This escape fraction is an order of 



magnitude large r than what has been measured in the 
Galaxy (~6 %: iBland-Hawthorn fe Malonevlll999D and 
dwar fs to date (~3 %: iZastrow et al. 1 1201 lUBareer et al.l 
120121) . In passing, we note that 3D simulations of UV ra- 
diative transfer in superwinds indicate that highe r values 
are possible in extreme cases (jYajima et al.l[2009| ). Thus, 
it is quite unlikely that the M82 cap clouds are ionized 
by photons from the M82 starburst regions. We consider 
shock ionization to be a more likely explanation, as we 
discuss in the next section. 

4.2. Shock ionization 

In order to expla in the Ha emission in the cap, 
ILehnert et ail (|1999f ) first suggested that the M82 su- 
perwind can drive a shock into the underlying gas cloud. 
The gas metallicity, which is presently unknown, is likely 
to fall in the range 0.1-1 Zq depending on whether the 
gas is infalling or entrained with the wind flow. The 
threshold metallicity for infalling gas at the present epoch 
appears to be close to O. IZ^ in all observations of the 
Local Universe to date ([Ryan- Weber et all [2009). We 
regard this value as a lower limit because the outer 
H I envelope in the M81 group appears to be material 
stripped from the outer disk of one or more galaxies 
(|Chvnoweth fc Langstonl 120071 ). The cap is unlikely to 
be entrained gas from the disk because dense material is 
broken up very quic kly by Rayleigh- Taylor instabilities 
(|Cooper et al.l I2009T) . as observed in the wind filaments 
close to the disk. 

We compare the observed line ratios with those cal- 
culated from a fas t shock model (v B i > 200 km s _1 ) of 
lAllen et al.l (|2008D . The [N II] /Ha ratios in the shock 
model are larger than 0.3, while the observed values are 
mostly less than 0.3 (Figure [3] (b)). But this may reflect 
the lower expected metallicity in the cap. Fast shocks in 
a low metallicity gas may be able to explain the observed 
[N II] /Ha ratios. The computed [N II] /Ha ratios fall in 
the range 0.1 to 0.3, which is similar to the observed 
values. But if fast shock excitation is dominant, the ob- 
served [N II] /Ha ratios should correlate with [S II] /Ha, 
because both ratios increase in lock step with an increase 
in shock velocity (Figure [3] (c)); however, no such corre- 
lation between these ratios is found. A more compelling 
argument against fast shocks is the kinematically 'cold' 
line emission observed across the 'cap' region, an issue 
we return to below. Thus fast shocks are unlikely to be 
the dominant ionization source of the cap today. 

Next we compare the observed ratios with those calcu- 
lated from slow to intermediate s hock models, i.e.. 40 km 
s" 1 <v sl < 130 km s" 1 . given bv lShull fc McKeel dl97l) 
(Figure [3] (a)). In this shock velocity range, [N Il]/Ha ra- 
tios increase as shock velocity increases, while [S II] /Ha 
ratios neither increase nor decrease monotonically. In 
Figure [3] (a), the observed points at the cap knots are 
distributed along the line of slow shock model. The cal- 
culated [N II] /Ha and [S II] /Ha ratios are 0.10 and 0.58 
at u Sj i = 60 km s^ 1 , and 0.27 and 0.44 at v sA = 80 
km s _1 , respectively. Therefore, shock velocities of 60 
km s _1 and 80 km s _1 can also explain the observed 
[N II] /Ha and [S II] /Ha ratios. Model comparisons for 
shock velocities higher than 80 km s _1 are ruled out. 

The superwind is powered by kinetic energy from some 
combination of stellar winds, radiation pressure and su- 
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Figure 5. Radial pressure profile along the M82 minor axis. 
Crosses, open squares, and an open circle r epresent the thermal 
pressure of ionized gas from Hcckman et al. (1990), Yoshida ct al. 
I|2011f ). and this study, respectively. The temperature of ionized 
gas is assumed to be 10 4 K. A filled squar e and a filled circle rep- 
resent the thermal pressure of h ot gas from Strickland & H eckmanl 
H2007I) and ILehnert et al] H999), respectively. The filling factor of 
hot gas is assumed to be unity. Thermal pressures are calculated 
by 2nkgT in both cases, where n and are the electron density 
of ionized gas and Boltzmann's constant, respectively. The slope 
of the dashed line is given by r~ 2 , where r is the distance from the 
M82 center. 

pernovae in the M82 starburst region. Here we ex- 
amine whether the momentum is enough to produce 
the observed shock. We compare the inferred range 
of shock velocities (40-80 km s -1 ) with the superwind 
velocity of a sphe r ically- symmetric model for M82 by 
iChevaiier fc Cleggl (j!985[) . Their model fits well with 
the observed thermal pressure profile of M82 within ~1 
kpc, and the pressure profiles at the larger radii than 
0.5 kpc are consist ent with the r~ 2 depen dence expected 
for ram pressure (Hcckma n et al] IT990). Using their 
model, the gas density, wind velocity, and ram pres- 
sure at the distance of M82 cap are estimated as 4 x 
10~ 6 cm" 3 , 5600 km s -1 , and 2.1 x 10~ 12 dyne cm" 2 , 
respectively. Since the observed electron density is 1.0 
cm~ 3 at the knot H-C, the thermal pressure of ionized 
gas is calculated as Ink^T — 2.8 x 10~ 12 dyne cm~ 2 
with the assumption of T — 10 4 K. This indicates that 
the model ram pressure and the observed thermal pres- 
sure are well balanced, and the observed pressure fol- 
lows r~ 2 law even at the distance of the M82 cap, 11.6 
kpc (Figure [SJ. Furthermore, given that some of the 
shock motion is expected t o generate local turbulent mo - 
tions of the same order ()Bland- Hawthorn et al.l I2007D . 
the measured narrow emission line widths (< 100 km 
s^ 1 FWHM, iDevine fc Ballvl I1999D are consistent with 
a slow shock. Therefore, it is quite likely that the cap 
clouds are ionized by a slow shock produced by the M82 
superwind. 

Future deep observations of the M82 cap, e.g., 
[O I]A6300, will enable us to directly confirm its dom- 
inant ionization source. In the case of a slow shock, 
the [O Il/Ha ratio l arger than 0.1 is to be expected 
(jShull fc McKedll979D . while a smaller rat io is evidence 
for p hotoionization by massive stars (e.g., IKewlev et al] 
12001 . 

4.3. Hot gas heating 

Diffuse X-ray emiss i on was detected at t he M82 cap 
(|Devine fc Ballvl 119991: ILehnert et all 119991: iTsuru et all 



2007). The electron density, temperature, and ther- 
mal pressure of hot gas are 5.4 x lO" 3 /^ 1 ^ 2 cm~ 3 , 9 

xlO 6 K, and 1.3 x 10~ n .f x ^ 2 dyne cm~ 2 , respectively, 
where fx is the filling factor of hot gas. However, a 
slow shock is not the main heating mechanism, because 
these can only heat to ~ 10 5 K. Additionally, thermal 
pressures of ionized and hot gas are inconsistent at the 
M8 2 cap, though they are consistent at the M82 cen - 
ter (jHeckman et al.lll990t [Strickland fc Heckmanl 120071) . 
Hence, another heating source for hot gas is required. 

A reasonable explanation for the hot gas is found when 
considering a 2-component shock model. An H I cloud 
in isolation develops a core-halo structure w here the core 
is dense and the hal o is relatively diffuse (Field 1965; 
ISternberg et al.ll2002D . In this interpretation, the X-ray 
emission is produced by the superwind triggering a shock 
in the diffuse gas su rrounding the knots in the cap (see 
ILehnert et all 1 19991 ). To produce gas with T = 9 x 10 6 
K requires a fast shock whose shock velocity v St 2 is 820 
km s _1 . But the 2-component shock by the same super- 
wind cannot explain the discrepancy between the ther- 
mal pressures of these gas phases. 

A 2-component shock by different superwinds from the 
M82 starburst region may be the case for the M82 cap. 
In this model, the slow shock by the present superwind 
ionizes the cap clouds, while hot gas observed today 
was produced by a fast shock driven by a past super- 
wind outburst. What is important is that the cooling 
timescales of these gas phases are drastically different: 
w 10 4 yr and ~ 10 8 yr for ionized and hot gas, respec - 
tively (jOsterbrock fc Ferlandll200l ILehnert et al.lll999D . 
The history of this model is as follows: < 10 8 yrs ago, 
a fast shock driven by a past superwind produced hot 
ionized gas; after « 10 4 yrs, the ionized gas cooled while 
the hot gas continued to emit X-rays, and now the slow 
shock ionizes the cap clouds again to produce the Ha 
and X-ray emission from the cap. 

5. CONCLUDING REMARKS 

It is undoubtedly true that the physics of super- 
winds is complicated. However, we are able to deduce 
some basic properties of the w ind. As discussed by 
ISharp fc "B land-Ha wthorn! (|2010D . there is a well deter- 
mined series of events that lead to a superwind taking 
hold. We believe that the new observations support the 
developing paradigm. 

When a critical surface density of dense molecular 
clouds is reached in the galaxy, massive stars are born. 
These stars evolve rapidly and their strong UV radiation 
fields produce a warm gaseous medium that encircles the 
remaining molecular clouds. After a few million years, 
the cores of the most massive stars collapse leading to 
multiple supernova explosions. This huge impulse of me- 
chanical energy heats the diffuse medium to extremely 
high temperatures (T ~ 10 8 K) where it expands to form 
a powerful superwind. (If the mechanical energy or UV 
radiative energy were entirely absorbed by the dense gas, 
most of the energy would simply be re-radiated as IR 
emission.) The hot flowing gas entrains cooler gas from 
the disk with the flow; we know this because of well de- 
fined rotation of t he entrained filaments about the wind 
axis (|Grevd 12004V The entrained gas is very clumpy 
because it is mostly entrained from the surviving dense 
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clouds in the disk ()Cooper et al.l 120081 ) . 

M82 is engulfed by a large H I cloud complex and some 
of this material appears to be accreting onto the dwar f 
close to the minor axis (jChynoweth fc Langstonl (2007 h 
The cap material is almost certainly supplied by infalling 
gas. The cap ionization almost certainly arises from the 
expanding superwind interacting with infalling gas. The 
spatial coincidence of the Ha and X-ray emission may be 
explained with the model of a 2-component shock by dif- 
ferent M82 superwinds. Further constraints on the cap 
may come from UV absorption line spectroscopy using 
distant quasars as a background light source. The ex- 
pected warm-hot medium should be visible in low to in- 
termediate ionization states of C, N, O, Ne and Fe, even 
for g as in a non-equilibrium state (|Gnat fc Sternberg] 
2004). Our new highly resolved observations show that 
the cap is very clumpy presumably as a consequence of 
the wind-cloud interaction. Our expectation is that the 
cap will be disrupted in a shock crossing time of about 
~ 1 Myr. 

The fact that the outer cap is presently ionized by lo- 
cal shocks, and not by nuclear UV radiat i on, su pports 
the findings of ISharp fc Bla nd-Hawthornl (|2010f ). Un- 
like AGN-driven wind filaments which are all found to 
be ionized by nuclear UV radiation, the starburst wind 
filaments are ionized by local shocks far from the nu- 
cleus. This is easily understood in terms of a 'starburst' 
because the hot young stars must evolve to supernovae 
before the wind gets going, and therefore few remain to 
ionize gas clouds in the direction of the flow. We fully an- 
ticipate that future deep optical, UV and X-ray imaging 
and spectroscopy will reveal further details about the cap 
region, and in turn about the nature of the superwind in 
M82. 

This work is based on data collected at Subaru Tele- 
scope, which is operated by the National Astronomical 
Observatory of Japan. We thank M. Yoshida for useful 
discussions. This work was financially supported in part 
by the Japan Society for the Promotion of Science (Nos. 
17253001, 19340046, 23244031, and 23654068). 
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